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Abstract 

Tsetse flies use olfactory and gustatory responses, through odorant and gustatory receptors (ORs and GRs), to interact with 
their environment. Glossina morsitans morsitans genome ORs and GRs were annotated using homologs of these genes in 
Drosophila melanogaster and an ab initio approach based on OR and GR specific motifs in G. m. morsitans gene models 
coupled to gene ontology (GO). Phylogenetic relationships among the ORs or GRs and the homologs were determined 
using Maximum Likelihood estimates. Relative expression levels among the G. m. morsitans ORs or GRs were established 
using RNA-seq data derived from adult female fly. Overall, 46 and 14 putative G. m. morsitans ORs and GRs respectively were 
recovered. These were reduced by 12 and 59 ORs and GRs respectively compared to D. melanogaster. Six of the ORs were 
homologous to a single D. melanogaster OR (DmOr67d) associated with mating deterrence in females. Sweet taste GRs, 
present in all the other Diptera, were not recovered in G. m. morsitans. The GRs associated with detection of CO2 were 
conserved in G. m. morsitans relative to D. melanogaster. RNA-sequence data analysis revealed expression of GmmORlS 
locus represented over 90% of expression profiles for the ORs. The G. m. morsitans ORs or GRs were phylogenetically closer 
to those in D. melanogaster than to other insects assessed. We found the chemoreceptor repertoire in G. m. morsitans 
smaller than other Diptera, and we postulate that this may be related to the restricted diet of blood-meal for both sexes of 
tsetse flies. However, the clade of some specific receptors has been expanded, indicative of their potential importance in 
chemoreception in the tsetse. 
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Introduction 

Trypanosomiasis management has been a longstanding devel- 
opment preoccupation in sub-Saharan Africa, with tsetse fly 
control constituting the cornerstone in this efibrt [1]. Since aU 
tsetse species are able to transmit trypanosomes, the critical 
determinant of transmission is their obligate blood feeding. Tsetse 
flies select their hosts through visual and olfactory signals, a 
process that is mediated by olfactory and gustatory receptors. 
Tsetse flies navigate their environment by detecting and respond- 
ing to volatUes and non-volatile cues (odors and tastants). Artificial 
bait technologies, based on tsetse olfactory responses to natural 
cues and blends of synthetic versions that mimic those of their 
natural hosts in the field, have successfully been applied in tsetse 
control because of their relatively high specificity, low cost, 
community acceptability, and ability to slow down tsetse re- 
invasion from adjacent areas [2,3]. These technologies are 
environment friendly [4] , and applicable for riverine and savanna 
species of tsetse flies [5,6]. The attractants include various phenolic 
derivatives [7-9], carbon dioxide, acetone, l-octen-3-ol, and 



vertebrate host breath, skin and urine extracts [10-12]. Interest- 
ingly, l-octen-3-ol is a constituent of the chemical profile from 
Lantana camara, an invasive plant to which tsetse flies are attracted 
[13]. The response to olfactory cues has also been exploited in 
design of tsetse repellents [14,15]. The repellents include guaiacol 
(methylphenols), 5-octalactone and methylketones [16-18] and 2- 
methoxy-4-methylphenol [14]. Natural differential responses 
among tsetse species and even between sexes and allopatric 
populations of the same species have been observed [18-22], 
which have stimulated research and design to enhance the 
efficiencies of the existing attractant-based bait technologies, to 
develop new ones based on repellent blends ('push' tactics) from 
refractory animals, and to integrate these into 'push-pull' 
strategies. Different Glossina species exhibit different olfactory 
uniqueness' and this may partly account for the observed 
graduation of preferences for particular hosts. For instance, 
riverine tsetse species (such as G. Juscipes Jiiscipes, G. palpalis and 
G. tachinoides) prefer feeding on reptilian hosts compared to their 
savanna relatives (G. morsitans morsitans, G. pallidipes) that feed 
largely on ungulates and other large mammals [6]. Larvipostion 
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Author Summary 

Tsetse flies navigate their environments using cliemosen- 
sory receptors, wliicli permit them to locate hosts, mating 
partners, and resting and larviposition sites. The genome 
of G. m. morsitans was interrogated for coding genes of 
odorant receptors (ORs) and gustatory receptors (GRs) that 
express in antennae and maxillary palp, and detect the 
volatile and soluble chemical signals. The signals are then 
transmitted to the central nervous system and translated 
to phenotypes. Majority of these genes in G. m. morsitans 
were spread across different scaffolds, but a few were 
found to occur in clusters, which suggested possible co- 
regulation of their expression. The number of ORs and GRs 
were much reduced in the G. m. morsitans genome, 
including the apparent loss of receptors for sugar when 
compared to selected Diptera. There was also an apparent 
numerical expansion of some receptors, presumably to 
maximize on their restricted blood-meal diet. The annota- 
tion of the chemoreceptor package of G. m. morsitans 
provides a resource for investigating key activities of tsetse 
flies that could be exploited for their control. 

pheromones (n-pentadecane and n-dodecane) from exudates of 
mature larvae are also known to attract and induce gravid G. m. 
morsitans and G. m. centralis females to aggregate and deposit larvae 
[23]. Research on response to tastants in tsetse flies are limited, but 
point to their potential application in tsetse control [10,24]. In all, 
responses to odors and tastants in tsetse have established utility in 
tsetse control that can be augmented with better understanding of 
the molecular factors that underpin these responses. 

Molecular factors mediating the olfactory and gustatory 
responses in the tsetse flies are poorly understood. However, 
research on other insects indicates that the odors and tastants in 
the environment are generally detected in peripheral sensory 
neurons by distinct odorant and gustatory receptors (ORs and 
GRs) [25-28]. These receptors are divergent members of a 
superfamily characterized by seven transmembrane domains, and 
share low sequence conservation among them except at the C- 
terminus region that coincides with the seventh trans-membrane 
domain [29]. The ORs and GRs are thought to have evolved as 
parallel chemoreceptors across diverse organisms [26]. Each OR is 
expressed in olfactory receptor neurons (ORNs) within maxiUary 
palpi and antennae [25,30-32]. The ORs generally have multiple 
introns and are very divergent with poor structural conservation 
within and across insect orders and species [33,34], which 
potentially reflect diverse olfaction related preferences within the 
orders and species. However, a canonical co-receptor commonly 
referred to as Oreo remains highly conserved across insect orders 
[35-38], a phenomenon that may be associated with its role in 
proper tuning of odor specificity and activation necessary for 
appropriate signal transduction in the neurons [39] . The GRs are 
generally expressed in gustatory receptor neurons (GRNs) within 
gustatory organs [40] in response to soluble taste and contact 
pheromones [41,42]. However, some GRs are expressed in 
antennal dendrites and respond to carbon dioxide, potentially 
implicating them in olfaction [40,43]. The GRs are more 
conserved in sequence and structure than the ORs [44,45] 
probably due to comparatively smaUer search space among cues 
associated with GRs than ORs. The diversity among the ORs and 
GRs in tsetse can potentially shed light on the natural differential 
responses observed among them [12,17,18,20-29], with potential 
application in tsetse control. To improve or develop new 
approaches of vector management, an understanding of the 



molecular attributes of GRs and ORs and their potential roles in 
tsetse ecology is essential. 

This study was initiated to (1) comparatively annotate and 
catalogue ORs and GRs in G. m. morsitans (GMOYl.l), (2) 
establish evolutionary distance between G. m. morsitans ORs or 
GRs and those in especially D. melanogaster, and (3) examine 
relative expression of the ORs and GRs in the G. m. morsitans. The 
assembly has been estimated to be over 99% complete based on 
the software Core Eukaryotic Genes Mapping Approach 
(CEGMA) [46] and manually sequenced BACs data. The 
assembly is currently undergoing genome-wide manual curation 
and annotation by the International Glossina Genome Initiative 
(IGGI) consortium. 

Methods 

Retrieval and annotation of G. m. morsitans OR and GR 
gene models 

Coding sequences (CDS) of ORs and GRs in Drosophila 
melanogasta' were obtained from FlyBase5.13 [47] and used to 
isolate their respective homologs in the G. m. morsitans genome 
(GMOYl.l) at VectorBase [48] using tBLASTx algorithm [49]. 
Scaffolds encoding the homologs were searched for and retrieved 
at a cut-off e-value <1.0e-05. Whole transcriptome Alumina 84 
million RNA sequence reads generated from female G. m. morsitans 
[50] were mapped onto the scaffolds using default settings in CLC 
Genomics workbench suite Version 4.8 (CLC Bio, Aarhus, 
Denmark). Gene loci of putative Glossina homologs were curated 
in the scaffold sequences flanking the tBLASTx hits, and intron/ 
exons modeled using the RNA-seq mappings. The predicted gene 
models were viewed and edited using Artemis vl3.2.12 [51] 
where, intron/ exon boundaries were edited using motifs GT for 5' 
donor site, and AG for 3' acceptor site. The start codon (ATG) for 
each gene model was fixed at the 5' end and the reading frame 
terminated at the first of any of the stop codons (TAA, TGA, or 
TAG). Sequences shorter than average size of known insect ORs 
(370 aa) were marked as incomplete if they lacked start or stop 
codons. Sequences with poorly conserved functional domains were 
considered as pseudogenes. 

The homologs were validated through sequence-based searches 
for presence of ORs or GRs specific 7tm-6-olf-recpt or 7tm-7-olf- 
recpt [29,52] domains respectively. The homologs were probed for 
the domains using DELTA BLAST algorithm [53] against the 
conserved domains databases (CDD) [54], and presence of alpha 
helix trans-membrane domains validated using TMHMM server 
v2.0 [55]. Additionally, all the putative ORs or GRs were 
validated, using BLAST2GO analyses [56] against the non- 
redundant Swiss-Prot database [57]. The curated gene models 
were assigned annotation identifiers by comparing them with 
automated transcript feature models obtained from the Glossina 
community annotation portal at VectorBase [48] and edited using 
Artemis genome viewer tool [51]. The models without automated 
prediction matches and identifiers were manually built using the 
Artemis gene build tool window [5 1] and given unique temporary 
annotation identifiers. In this respect, features for gene, exons, 
mRNA, and CDS were created for such gene models. The 
Glossina gene models were assigned putative gene names where 
GmmOR* and GmmGR* were adopted for G. m. morsitans 
odorant receptors and gustatory receptors respectively (the asterisk 
(*) being an identifier number). The annotated gene model 
features were submitted to the VectorBase community annotation 
portal for G. m. morsitans [48] for integration into genome database; 
nevertheless, a list of annotated amino acid coding sequences is 
presented in supplementary Dataset S 1 , and a list of associated 
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Table 1. Annotated ORs and GRs in G. m. morsitans and other 
selected insect species. 



Insect 


ORs 


GRs 


Reference 


D. melanogaster 


60 (2)* 


60 03)* 


[25-27,29,42] 


G. m. morsitans 


46 (3) 


14 


This study. 


An. gambiae 


79 


76 


[52,76] 


Ae. aegypti 


100(31) 


79 


[74,77] 


Apis meliifera 


163 (11) 


10 (3) 


[41] 


Nasonia vitripennis 


225 (76) 


47 (11) 


[80] 


Camponotus floridanus 


352 (55) 


46 (17) 


[75] 


Harpegnathos sattator 


347 (30) 


17 (4) 


[75] 


Tribotium casteneum 


265 (76) 


220 (25) 


[78,79] 



Figures in parentheses are numbers incomplete genes and or pseudogenes of 
the receptors. 

*- in parentheses are alternatively spliced forms. 
doi:l 0.1 371/journal.pntd.0002663.t001 



gene identities in Table S2. The G. m. morsitans receptor repertoires 
were evaluated against those documented for D. melanogaster. 
Anopheles gambiae, Aedes aegfpti, Apis meliifera, Nasonia vitripennis, 
Camponotus Jloridanus, Harpegnathos saltator and Tribolium casteneum 
(references in Table 1). 

Phylogenetlc analyses of ORs and GRs in G. m. morsitans 
and selected Diptera 

Multiple Sequence Comparison by Log-Expectation (MUS- 
CLE) tool [58] was used to align GmmORs and GmmGRs with 
homologs in D. melanogaster, and the alignments edited using 
Jalview web-server [59]. The secondary structures in the 
alignments were predicted using JPred program [60]. Phylogenetlc 
cluster inference was done using Maximum Likelihood approach 
with best fitting Wheelan and Goldman+Freq (WAG+F) model 
[61], which was chosen as the best ranked from a panel of all 
amino acid model tests run in MEGA5 [62] . The initial tree was 
automatically generated and bootstrapped with 500 iterations. 
The evolutionary rate difference among sites was modeled using a 
discrete Gamma distribution (5 categories (+G, parame- 
ter = 4.2651)). The rate variation model allowed for some sites to 
be evolutionarHy invariable ([+1], 0.8705% sites). All positions with 
less than 95% site coverage were eliminated and branch nodes 
determination set at very strong. Evolutionary analyses were 
conducted using the MEGA5 suite [62] . 

Comparative analyses of expression profiles of G. m. 
morsitans ORs and GRs 

The expression profiles of G. m. morsitans ORs and GRs gene 
loci were determined using whole transcriptome 84 million 
Ulumina RNA-sequence reads [50]. The RNA-seq reads were 
mapped onto the G. m. morsitans ORs or GRs nucleotide coding 
sequences (CDS) in CLC Genomics Workbench (CLC Bio, 
Aarhus, Denmark) via RNA-seq analysis pipeline with default 
settings. The expression profiles were presented as reads per 
kilobase of exon model per million mapped reads (RPKM) for 
each receptor sequence [63]. 

Results 

Most of the gene loci of G. m. morsitans ORs and GRs were 
scattered amongst the scaffolds. Fifty percent of G. m. morsitans 



OR genes were encoded as single-copies on their respective 
scaffolds. The remainder were encoded in pairs or triplets per 
scaffold. Five G. m. morsitans OR loci (GmmOR6/7/8, 
GmmOR18/19, GmmOR22/25, GmmOR27/28 and 
GmmOR41/42) were located in tandem on their respective 
scaffolds. Similarly, five G. m. morsitans GR genes clustered on a 
single scaffold. The rest were encoded as single-copies on their 
respective scaffolds. All G. m. morsitans GR loci were annotated as 
complete genes. 

Gene models for G. m. morsitans OR and GR and their 
annotation 

Numbers of OR and GR gene loci recovered in G. m. morsitans, 
relative to those published in other insects are summarized in 
Table 1 . Similar to most insects, the G. m. morsitans has more ORs 
loci than GRs loci, with the exception of D. melanogaster where the 
numbers are equal. However, the G. m. morsitans ORs are fewer 
than those documented in all the insects evaluated, including D. 
melanogaster. A similar trend was exhibited in G. m. morsitans GRs, 
except in relation to A. meliifera. Annotation of G. m. morsitans ORs 
and GRs are summarized in Table 2. The lengths of G. m. morsitans 
ORs varied between 260 and 541 amino acids, while those of G. m. 
morsitans GRs ranged from 309 to 5 14 amino acids. The number of 
exons ranged between two and eight or 12 in GRs and ORs 
respectively. The predicted genome structures are given in Figure 
SI. The frequency of detectable trans-membrane domains was 
also variable, with proteins having six trans-membrane domains 
representing about one half of all genes. The G. m. morsitans ORs 
(57%, 26 out of 46) were homologous to nine D. melanogaster ORs. 
Similarly, most of the G. m morsitans GRs (57%, 8 out of 14) were 
homologous to three D. melanogaster GRs genes. The remainder of 
the G. m. morsitans GRs had one-to-one homology with a single D. 
melanogaster specific homolog. Reciprocal blasts onto non-redun- 
dant protein databases for both G. m. morsitans ORs and GRs are 
summarized in Supplementary material - Table SI). GmmGR3 
and GmmGR4 were also homologous to An. gambiae orthologs, 
while GmmGR5, GmmGR8 and GmmGR13 had homologs to 
genes in other Drosophila species. The G. m. morsitans ORs 
pseudogenes were scanty, representing 7% of the ORs genes 
recovered. Only GmmOR5 had alternative splice variants. The 
7tm-6-olfct-rcpt domain was detected in aU G. m. morsitans ORs, 
and the 7tm-7-chem-rcpt domain was detected in five ORs 
(GmmOR17, GmmOR21, GmmOR24, GmmOR38 and 
GmmOR39). The 7tm-7-chem-rcpt domain was also detected in 
all the G. m. morsitans GRs. 

Phylogenetlc analysis of G. m. morsitans ORs and GRs 
with other insects 

Phylogenetlc relationships between G. m. morsitans ORs and GRs 
and their counterparts in D. melanogaster are summarized in 
Figure 1 . Most of the G. m. morsitans ORs and GRs clustered with 
their respective ORs and GRs orthologs with a bootstrap support 
of over 80%. The G. m. morsitans OR 14, OR 15 and OR 16 were 
homologous to a drosophila larvae receptor, Or45a. The G. m. 
morsitans co-receptor (Oreo) (GmmORl) had 100% bootstrap 
support homology to D. melanogaster homolog, Or63b, and was a 
single copy in the genome, similar to other insects investigated 
(data not shown). There was an expanded cluster of ORs in G. m. 
morsitans (GmmOR41-46), relative to a single D. melanogaster 
homolog, Or67d (Figure lA), which also had multiple copies in 
An. gambiae, Cu. quinquefasciatus, Ae. aegypti, Tribolium casteneum (Data 
not shown). The G. m. morsitans and D. melanogaster GRs clustered 
into four groups (Figure IB). Four G. m. morsitans GRs (GmmGRl- 
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Table 2. Annotations of odorant and gustatory receptor genes in G m. morsitans and their homologs in D. melanogaster. 





G. m. morsitans Genes 


Length (AA) 


Exons 


TMMs 


Gene ID 


Dmel orthologs/Accesslon Number 


GmmOR! 


521 


8 


6 


GMOY005610 


DmOr83b/CG 10609 


GmmOR2 


394 


3 


7 


GMOY005796 


DmOr2a/CG3206 


GmmOR3 


387 


3 


3 


GMOY004772 


DmOr19a/CG18859 


GmmOR4 


384 


2 


7 


TMP_0r4 


DmOr59a/CG9820 


GmmORS* 


442 


4 


5 


GMOY012018 


DmOr33b/CG16961 


GmmOR6 


387 


4 


5 


GMOY009475 


DmOr42b/CG 12754 


GmmOR7 


406 


3 


6 


TMP_0r7 


DmOr42b/CG 12754 


GmmORS 


389 


4 


6 


TMP_0r8 


DmOr42b/CG 12754 


GmmOR9 


409 


3 


6 


TMP_0r9 


DmOr42b/CG 12754 


GmmORlO 


444 


3 


6 


TMP_Or10 


DmOr46a/CG33478 


GmmORll 


341 


3 


6 


GMOY010761 


DmOr46a/CG33478 


GmmOR12 


340 


3 


3 


GMOY009271 


DmOr94b/CG 17241 


GmmOR13 


391 


6 


6 


GMOY003312 


DmOr82a/CG31519 


GmmORM 


341 


3 


6 


GMOY001365 


DmOr45a/CG1978 


GmmORlS 


446 


4 


7 


TMP_0rl5 


DmOr45a/CG1978 


GmmORie 


387 


4 


6 


TMP_0rl6 


DmOr45a/CG1978 


GmmOR17 


541 


12 


8 


GMOY005386 


DmOr69a/CG33264 


GmmORlS 


420 


8 


6 


TMP_0rl8 


DmOr63a/CG9969 


GmmOR19 


385 


8 


7 


GMOYO 12322 


DmOr63a/CG9969 


GmmOR20# 


269 


7 


6 


TMP_Or20 


DmOr85b/CG11735 


GmmOR21 


465 


5 


2 


GMOYO 11399 


DmOr83a/CG10612 


GmmOR22# 


296 


4 


5 


TMP_Or22 


DmOr49a/CG13158 


GmmOR23 


331 


4 


5 


TMP_Or23 


DmOr85b/CG11735 


GmmOR24 


388 


3 


6 


GMOYO 10839 


DmOr85c/CG17911 


GmmOR25 


385 


3 


6 


GMOY012357 


DmOr56a/CG12501 


GmmOR26 


418 


4 


5 


TMP_Or26 


DmOr85b/CG11735 


GmmOR27 


415 


3 


6 


GMOYCC8038 


DmOr67c/CG14156 


GmmOR28# 


260 


2 


7 


TMP_Or28 


DmOr92a/CG17916 


GmmOR29 


438 


3 


4 


TMP_Or29 


DmOr67a/CG12526 


GmmORBO 


361 


3 


6 


TMP_Or30 


DmOr67a/CG12526 


GmmOR31 


435 


7 


5 


TMP_0r31 


DmOr24a/CGn767 


GmmOR32 


450 


5 


7 


GMOY005084 


DmOr13a/CG12697 


GmmOR33 


353 


6 


5 


GMOY005479 


DmOr49b/CG 17584 


GmmOR34 


360 


7 


4 


GMOYO 11902 


DmOr30a/CG13106 


GmmOR35 


392 


5 


6 


TMP_Or35 


DmOr43a/CG1854 


GmmOR36 


343 


7 


6 


TMP_Or36 


DmOr43a/CG1854 


GmmOR37 


430 


4 


4 


TMP_Or37 


DmOr74a/CG13726 


GmmOR38 


371 


5 


6 


TMP_Or38 


DmOr47b/CG 13206 


GmmOR39 


403 


3 


6 


GMOY004392 


DmOr88a/CG14360 


GmmOR40 


284 


5 


6 


GMOY012356 


DmOr56a/CG12501 


GmmOR41 


386 


4 


6 


GMOY006480 


DmOr67d/CG14157 


GmmOR42 


386 


4 


5 


GMOY006479 


DmOr67d/CG14157 


GmmOR43 


389 


4 


5 


TMP_Or43 


DmOr67d/CG14157 


GmmOR44 


390 


4 


6 


GMOY006265 


DmOr67d/CG14157 


GmmOR45 


385 


4 


7 


GMOY007896 


DmOr67d/CG14157 


GmmOR46 


348 


4 


3 


GMOY003305 


DmOr67d/CG14157 


GmmGRl 


425 


3 


6 


GMOY007472 


DmGr21a/CG13948 


GmmGR2 


514 


7 


6 


GMOY011510 


DmGr22b/CG31931 


GmmGR3 


425 


6 


6 


TMP_Gr5 


DmGr21a/CG13948 
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Table 2. Cont. 



G. m. morsitans Genes 


Length (AA) 


Exons 


TMMs 


Gene ID 


Dmel orthologs/Accesslon Number 


Ol 1 1[ 1 lUn^ 


496 


g 


g 


UIVIW I UUOUU 1 






467 


J 


7 


UlVIU I uu^zu/ 




GmmGR6 


443 


4 


g 


UlVIU I U 1 1 O 1 J 


UmUl ZdD/lwL] 1 D/oo 


GmmGR7 


402 


3 


7 


GMOY006209 


DnnGr28b/CG13788 


GmmGRS 


407 


2 


6 


TMP_Gr4 


DmGr22e/CG31936 


GmmGR9 


348 


5 


4 


GMOYC11903 


DmGr2a/CG18531 


GmmGRlO 


458 


4 


7 


GMOY003231 


DmGr33a/CG17213 


GmmGRll 


450 


3 


6 


TMP_Gr3 


DnnGr22b/CG31931 


GmmGR12 


375 


2 


8 


TMP_Gr2 


DmGr32a/CG14916 


GmmGR13 


457 


2 


6 


TMP_Grl 


DmGr22b/CG31931 


GmmGRM 


309 


3 


6 


TMP_Gr6 


DmGr22b/CG31931 



GmmOR - Glossina morsitans morsitans ordorant receptor; GmmGR- G. m. morsitans gustatory receptor; TMM- Trans-membrane helices; GMOY - Glossina morsitans Yale 
strain; TMP_Or - Provisional odorant receptor ID; TMP_Gr - Provisional gustatory receptor ID; DmOr- Drosophila melanogaster odorant receptor; DmGR- D. melanogaster 
gustatory receptor; 

*- longest alternative splice variant in locus 0R5; 
^- pseudogene. 

doi:1 0.1 371 /journal.pntd.0002663.t002 



4) clustered with homologs of CO2 receptors, Gr21a and Gr63a in 
D. melanogaster, GmmGR6-7 and GmmGR 14, though distantly, 
clustered with an unusual splice variant DmelGr28a/28b; 



GmmGRS, 8-12 were homologous to bitter taste-related sensors 
in D. melanogaster, and GmmGR 13 clustered distantly to Dmel- 
Gr58a/58b homologs, whose functions are unknown. 




Figure 1. Phylogenetic analyses of ORs or GRs in G. m. morsitans and selected Diptera. (A) Maximum likelihood (ML) tree for GmmORs and 
DmelOrs; branches annotated blue is an expanded clade orthologous to DmelOr67d; purple branches is the clade orthologous to DmelOr45a; and 
green branches indicate the orco cluster. (B) Maximum likelihood tree for GmmGRs and DmelGRs. In both trees, blue labels are D. melanogaster 
receptors and red labels G. m. morsitans receptors (green labels are An. gambiae C02 receptors). Phylogenetic cluster inferences were deduced using 
Maximum Likelihood approach with best fitting Wheelan And Goldman+Freq (WAG+F) model [59]. Evolutionary analyses were conducted using 
MEGA5 suite [60]. 

doi:1 0.1 371/journal.pntd.0002663.g001 
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Figure 2. Glossina chemoreceptor expression abundances by RNA-seq data in RPKM. (A) Expression abundances of GmmORs. There was 
no sufficient data to support profiles for GmmORS, 1 1, 25, 31, and 39. GmmORI 5 had abundant transcriptome data of 90.746% relative to sequence 
reads that mapped onto GmmORs. [B] Expression abundances of GmmGRs. Expression profiles for GmmGR6, GmmGR1 1 and GmmGRI 3 were not 
detected; GmmGR2 and GmmGR3 accounted for 40% and 36% respectively of the total transcripts considered for GRs. 
doi:10.1371/journal.pntd.0002663.g002 



Relative expression profiles of G. m. morsitans ORs and 
GRs 

Relative expression profiles of the G. m. morsitans ORs and GRs 
gene loci are summarized in Figure 2. Among the G. m. morsitans 
ORs, expression of GmmORI 5 was surprisingly most predominant, 
accounting for more than 90% of the total RNA-sequence data 
supporting expression of the ORs. GmmORlS is homologous to 
Or45a gene in D. metanogaster. About 5% of RNA-sequence data 
provided supporting evidence for expression of GmmOR2, 
GmmORI (Oreo homolog), GmmOR43 and GmmORS. Expressions 
of GmmORS, GmmORI 1, GmmOR25, GmmORSl, and GmmOR39 
were not detected in the available RNA-sequence dataset 
(Figure 2A). Amongst the GRs, GmmGRI -4 had the best RNA- 
sequence data expression support (Figure 2B). 

Discussion 

Specific groups of the G. m. morsitans ORs and GRs were 
clustered within selected scaffolds. Similar clusters of genes 



performing common and related functions have been observed 
among chemosensory genes in D. melanogaster [41,42,44], and more 
recently among twelve G. m. morsitans major mUk proteins 
associated with lactation [50]. Since genes within clusters are 
generally co-regulated and can lead to joint gene expression 
[29,34,64], die individual clusters of ORs and GRs might be 
under common regulatory mechanisms and in response to 
common or related stimuli. The ORs and GRs in G. m. morsitans 
were fewer than those documented in most insects evaluated 
(Table 1) [65,66]. Additionally, specific ORs and GRs in D. 
melanogaster (nine and three ORs and GRs respectively) appear to 
have been expanded in G. m. morsitans, representing more than half 
of the chemoreceptors. 

The factors underlying the apparent reductions and expansions 
of these receptors in the tsetse are unknown. However, it can be 
postulated that the obligate blood feeding of the tsetse fly 
(restricted to vertebrate hosts) relative to D. melanogaster (with 
expansive fruit species hosts) might have necessitated evolutionary 
selection for specific chemoreceptor loci relevant to discriminate 
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among limited host choices. We know also that environmental 
factors can determine host choice, as tsetse have been shown to 
have an acquired preference to specific hosts encountered early in 
life [67]. Notably, other blood-feeders, such as mosquitoes also 
seek a variety of plant sources for sugar as energy source, while 
tsetse flies derive their energ)' from the amino acids proline and 
alanine [68]. The G. m. morsitans OR 15 (GmmORlS) accounted 
for more than 90% of the OR expression data. This OR is 
homologous to DmelOr45a, whose product has been, associated 
with an escape response in D. melanogaster larvae [69] . The function 
of this OR in tsetse was not determined; nonetheless it is notable 
that the source of RNA sequence data was a reproductively active 
adult female. Hence, it is possible that the GmmOR15 is in some 
way associated with larval activity. 

Similarly, the GmmGRl-4 cluster was most prominent among 
the GRs homologous to CO2 receptors in D. melanogaster. These GRs 
may be associated with host seeking and may have a duplicate role 
in olfaction. These receptors may putatively be associated with 
attractive responses elicited by the savanna tsetse species, including 
G. m. morsitans [1 0] . From the foregoing, it is evident that tsetse seems 
to prioritize and invest on a select few chemoreceptor genes towards 
their adaptive behaviors. Indeed, a heavy investment in specific 
genes is not uncommon in insects [70-73]. The G. m. morsitans ORl 
(homologous to Oreo) was the most conserved amongst the G. m. 
morsitans ORs, not surprising since such conservation has been 
observed in other insects [74] probably due to its critical role in 
modulating responses of the other receptors. 

In conclusion, when examined against other blood feeders, which 
also take sugar sources from plants (e.g. An. gamhiae and Ae. aegfipti), 
the G. m. morsitans hafi a reduced ri'p(;rt()ire of ORs and GRs genes. 
There is a complete loss of receptors for sugar, and a heavy 
investment in some chemoreceptors, such as those associated with 
detection of CO2. These observations offer opportunities to develop 
control tools exploiting these unique adaptations. 

Supporting Information 

Dataset SI Annotated amino acid sequences of Glossina 
morsitans morsitans ORs and GRs. Each sequence name is 

References 

1. Elliott M,Jancs Nl% Potter C (1978) The Future of Pyrethroids in Insect Control. 
Annual Review of Entomology 2.^: 443—469. 

2. Vale GA, Lovemore DF, Flint SCG (1988) Odour-baited targets to eontrol 
tsetse flies, Glossina spp. (Diptera: Glossinidae) in Zimbabwe. BuD Ent Res 78: 31- 
49. 

3. Mangwiro TNC, Torr SJ, CoxJR, HoUoway M (1999) The efficacy of various 
pyrethroid inseetieides for use on odour-baited targets to eontrol tsetse. Med Vet 
Ent 13: 31,S-323. 

4. Allsopp R (2001) Options for vector eontrol against trypanosomiasis in Africa. 
Trends Parasitol 17: 1,S-19. 

,'5. Gibson G, Forr SJ (1999) Visual and olfaetor)' responses of haeinatophagous 
Diptera to host stimuli. Med Vet Entomol 13: 2-23. 

6. Omolo MO, Hassanali A, Mpiana S, Esterhuizen J, Lindh J, ct al. (2009) 
Prospects for developing odour baits to control Glossina fuseipes spp., the 
major vector of human African trypanosomiasis. PLoS Ncgl Trop Dis 3: e4 
35. 

7. Den Otter CJ (1991) Olfactory responses of tsetse flies to phenols from buffalo 
urine. Physiol Entomol 16: 401-410. 

8. Den Otter CJ, Van der (iocs van Naters \VM (1993) Responses of individual 
olfactory cells of tsetse flies {Glossina m. morsitans) to phenols from cattle urine. 
Physiol Entomol 18: 43-49. 

9. Saini RK, Hassanah A, Andoke J, Ahuya P, Ouma WP (1996) Identification of 
major components of larviposition pheromone from larvae of tsetse flies Glossina 
morsitans morsitans Westwood and Glossina morsitans centralis Machado. J Chem 
Ecol 22: 1211-1220. 

10. BOGNER F (1992) Response properties of C02-sensitive receptors in tsetse flies 
(Diptera: G/ossina paipaiis). Physiol Entomol 17: 19-24. 

1 1. Vale GA (1980) Field studies of the response of tsetse flies ((}lossinidae) and other 
Diptera to carbon dioxide, acetone and other chemicals. Bull Entomol Res 70: 
563-570. 



followed by annotation identity (previous temporary identity), 
scaffold number and the coordinates within the scaffold. 

(PDF) 

Figure SI Glossina m. morsitans ORs and GRs genome 
structure. Most of the gene loci were encoded as singlets on their 
respective scaffolds. Some loci were encoded in tandem in their 

respective scaffolds suggestive of possible joint expression regula- 
tion. All genes had multiple exons ranging from two to eight or 1 2 
in GRs or ORs respectively. 

(PDF) 

Table SI Reciprocal blast of G. m. morsitans ORs and 
GRs to non-redundant protein database at NCBI. The 

names of annotated gene loci, scaffold identity, gene location 
within the scaffold, in bracket (*) refers to the coding strand where 
(— ), reverse and (-I-), forward strands. The reciprocal blast gives the 
distribution of homology to other insects including D. melanogaster. 
(PDF) 

Table S2 List of Glossina m. morsitans gene names and 
their associated identities. The vectorbase identities 
(GMOY ) has since replaced the Glossina community annota- 
tion identities (TMPC) . . . .) that were computationally derived. Both 
annotation identiti(-s and phylomedb identities can be used as seed 
queries to retrieve their related Glossinia phylomedb 182 trees 
(www.phylomedb.org/pq = user/28). 
(PDF) 

Acknowledgments 

Wc arc grafclul for the considerable technical support from the 
International Glossina Genome Initiative (IGGI) consortium, which also 
made the genome sequence available. 

Author Contributions 

Conceived and designed the experiments: GFOO POM .^C DKM. 
Performed the cxperimenls: GFOO. Analyzed the data: GFOO HMR 
SRtjN. Contributed reagents/maLerials/analysis tools: POM AC DKM. 
Wrote the paper: GFOO POM AC SRGN HMR DKM. Designed Glossina 
morsitans morsitans gene name/ symbol nomenclature: HMR. 



12. Vale GA, Hall DR (198,^)b) Fhr use of l-oeten-3-ol. aertonr and carbon dioxide 
to improve baits for tsetse flies, (Uossina spp. (Diptera: Glossinidae), to host odour. 
Bull Entomol Res 7.'"): 219 231. 

13. Syed Z, Guerin PM (2004) Tsetse flies are attracted to the invasive plant Lantana 
comma. Journal of insect physiology 50: 43—50. 

14. Saini RK, Hassanali A (2007) A 4-alkyl-substituted analogue of guaiacol shows 
greater repellency to savannah tsetse [Glossina spp.). J Chem Eeol logy 33: 985— 
995. 

15. Willemse L, Fakken VV (1994) Odor-indueed host location in tsetse flies (Diptera: 
Glossinidae), J Med Entomol 31: 77,')-791. 

16. Echane M (2005) lire Biology of Blood-Sircking in Insects. Second, (Cambridge, 
UK: CUP, 

17. Gikonyo NK, Hassanah A, Njagi PC, Gitu PM MJ (2002) Odor composiuon of 
preferred (buffalo and ox) and non-preferred (waterbuck) hosts of some Savanna 
tsetse flies. J Chem Ecol 28: 969-981. 

18. Gikonyo NK, Hassanali A, Njagi PGN, Saini RK (2003) Responses of Gtossina 
morsitans morsitans to blends of electroantennographically active compounds in the 
odors of its preferred (buffalo and ox) and non preferred (waterbuck) hosts. 
J Chem Eeol 29: 2331-2345. 

1 9. Mwangi MT, Gikonyo NK, Ndiege lO (2008) Repellent properties of delta 
octalactone against the tsetse fly Glossina mor.ntans morsitans.^ Insect Sci 8: 7-10. 

20. Vale GA (1985a) Flight as a factor in the host-finding behaviour of tsetse flies 
(Diptera: Glossinidae). Bull Entomol Res 70: 299-307. 

21. Vale GA, HaU DR, Gough AJE (1988) The olfactory responses of tsetse flies, 
Ghssina spp. (Diptera: Glossinidae), to phenols and urine in the field. Bull 
Entomol Res 78: 293-300. 

22. Mireji PO, Mabveni AM, Dube BN, Ogembo JG, Matoka CM, Mangawiro 
TNC (2003) Field responses of tsetse flies (Cilossinidae) and other Diptera to 
oils in formulations of deltamethrin. Int J Trop Insect Sci 23: 317- 
323. 



PLOS Neglected Tropical Diseases | www.plosntds.org 



7 



April 2014 | Volume 8 | Issue 4 | e2663 



Glossina morsitans morsitans Chemoreceptors 



23. Saini RK, Hassanali A, Aiidokc J, Ahuya P, Ouma \VP (1996) Identification of 
major components of larviposition phcromonc from larvae of tsetse flies Glossina 
morsitans morsitans Westwood and Glossina morsitans centralis Machado. J Chem Ecol 
22: 1211-1220. doi:10.1007/BF02266961. 

24. Den Otter CJ, Saini RK (1985) Pheromone perception in the tsetse fly. Entomol 
Exp Appl 39: 155-161. 

25. Clynr P[, Warr CO, Freeman MR, Lessing D, Kim J, Carlson JR, 1999 A novel 
familv ()[ di\ cri^( iit seven-membrane proteins; candidate odorant receptors in 
Drosoph'da. Neuron 22: 327-338. 

26. Vosshal LB,Amrein H, Morozow PS, Rzhetsky A, Axel A (1999) A spatial map 
of olfactory receptor expression in the Lhvsophila antenna. Cell 96: 725-736. 

27. Clyne PJ, Warr CG, Carlson JR (2000) Candidate taste receptors in Drosopkila. 
Science 287: 1830-1834. 

28. Ghaninia M, IgncU R, Hansson BS (2007) Functional classification and central 
nervous projections of olfactory receptor neurons housed in antennal trichoid 
sensilla of female vellow fe\'er mosquitoes, Aedes aegypti. EurJ Neurosci 26: 161 1- 
1623. 

29. Robertson HM, Warr CO, Carlson JR (2003) Molecular evolution of the insect 
chemoreecptor gene superfamily in Drosopkila melanogaster. Proc Nad Acad Sci 
USA 100: 14537-14542. 

30. Fuss W (2009) Does life originate from a single molecule? Chirality 2 1 : 299-304. 

31. Zwiebel LJ, Takken W (2004) Olfactory regulation of mosquito-host interac- 
dons. Insect Biochem Mol Biol 34: 645-652. 

32. Fuss SH, Ray A (2009) Mechanisms of odorant receptor gene choice in 
Drosophila and vertebrates. Mol Cell Neurosci 41: 101—112. 

33. de Bruyne M, Baker TC (2008) Odor detection in insects: volatile codes. J Chem 
Ecol 34: 882-897. 

34. Nozawa M, Nei M (2007) Evolutionary dynamics of olfactory receptor genes in 
Drosophila species. Proc Nad Acad Sci USA 104: 7122-7127. 

35. Krieger J, Raming K, Dewer YM, Bette S, Conzelmann S, et al. (2002) A 
divergent gene family encoding candidate olfactory receptors of the moth 
Helioihis virescens. EurJ Neurosci 16: 619-628. 

36. Jones WD, Nguyen TA, Kloss B., Lee KJ, Vosshall LB (2005) Functional 
conservation of an insect odorant receptor gene across 250 million years of 
evolution. Curr Biol 15: R119-R121. 

37. Vosshall TN, Vosshall LB (2009) Mechanisms in the Insect Olfactory System. 
Curr Opin Neurobiol 19: 284-292. doi:10.1016/j.conb.2009.07.015. 

38. Hansson BS, Stensmyr MC (2011) Review: Evolution of Insect Olfaction. 
Neuron 72: 698-711. 

39. PeUegrino M, Nakagawa T (2009) Smelling the difference: controversiEd ideas in 
insect olfaction. J Exp Biol 212: 1973-1979. 

40. Montell C (2010) A Taste of the Drosophila Gustatory Receptors. Curr Opin 
Neurobiol 19: 345-353. doi:10.1016/j.eonb.2009.07.001.A. 

41. Robertson HM, Wanner KW (2006) The ehemoreceptor superfamily in the 
honey bee. Apis mellifera: expansion of the odorant, but not gustatory, receptor 
family. Genome Res 16: 1395-1403. 

42. SAnchez-Gracia A, Vieira EG, Rozas J (2009) Molecular evolution of the major 
chemosensory gene families in insects. Heredity 103: 208-216. 

43. Fiala A (2007) Olfaction and olfactory learning in Drosophila: recent progress. 
Curr Opin Neurobiol 17: 720-726. doi:10.1016/j.conb.2007.1 1.009. 

44. McBride CS, Arguello JR, O'Meara BC (2007) Five Drosophila genomes reveal 
nonneutral evolution and the signature of host specialization in the 
chemoreecptor superfamily. (ienetics 177: 1395-1416. 

45. Gardiner A, Barker D, Butlin RK, Jordan WC, Ritchie MG (2008) Drosophila 
ehemoreceptor gene evolution: selection, specialization and genome size. Mol 
Ecol 17: 1648-1657. 

46. Parra G, Bradnam K, Korf I (2007) CEGMA: a pipeline to accurately annotate 
core genes in eukaryotic genomes. Bioinformatics (Oxford, England) 23: 1061— 
1067. 

47. MeQuilton P, Susan E, Pierre St, Thurmond J (2012) FlyBase Consortium 
FlyBase 101 - the basics of navigating FlyBase. Nucleic Acids Res 40: D706-14. 
doi:10.1093/nar/gkrl030. 

48. Megy K, Emrich SJ, Lawson D, Campbell D, Dialynas E, et al. (2012) 
VectorBase: improvements to a bioinformatics resource for invertebrate vector 
genomics. Nucleic Acids Res 40: 1-6. 

49. Altschul SF, Lipman DJ (1990) Protein database searches for multiple 
alignments. Proc Natl Acad Sci USA 87: 5509-5513. 

50. Benoit JB, Geoffrey MA, Veronika M, Tyler BK, Jana B, Qirui Zhang, et al. 
(2013) A novel highly divergent protein family from a viviparous insect identified 
by RNA-seq analysis: a potential target for tsetse tly-speeific abortifacients. 
PLCS Genetics (in Press). 

51. Carver T, Berriman M, Tivey A, Patel C, Bohme U, et al. (2008) Artemis and 
ACT: viewing, annotating and comparing sequences stored in a relationjil 
database. Bioinformatics 24: 2672-2676. 

52. Fox AN, Pitts RJ, Robertson HM, Carlson JR, Zwiebel LJ (2001) Candidate 
odorant receptors from the malaria vector mosquito Anopheles gambiae and 
evidence of down-regulation in response to blood feeding. Proc Nati Acad Sci 
USA 98: 14693-14697. 



53. Boratyn CM, Sehaffer A A, Agarwala R, Altschul SF, Lipman DJ, et al. (2012) 
Domain enhanced lookup time accelerated BLAST. Biology Direct 7: 12. 

54. Marchler-Bauer A LS, Anderson JB, Chitsaz F, Derbyshire MK, Deweese-Scott 
C, et al. (2011) CDD: a Conserved Domain Database for the functional 
annotation of proteins. Nucleic Acids Res 39: D225-229. 

55. Krogh A, Larsson B, von Heijne G, Sonnhaimner EL (2001) Predicting 
transmembrane protein topolo^gy with a hidden Markov model: application to 
complete genomes. J Mol Biol 305: 567-580. 

56. Conesa A, Gotz S, Ciarria-Clomez JM, I'erolJ, Talon M, et al. (2005) Blast2GO: 
a universal tool for annotation, visualization and analysis in functional genomics 
research. Bioinformatics (Oxlbrd, England) 21: 3674—3676. 

57. The UniProt Consortium (2013) Update on activities at the UniversEd Protein 
Resource UniProt in 2013 Nucleic Acids Res 41: D43-D47. 

58. Edgar Robert C. (2004), MUSCLE: multiple sequence alignment with high 
accuracy and high throughput. Nucleic Acids Res 32(5), 1792-97. 

59. Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ (2009) Jalview 
Version 2-a multiple sequence alignment editor and analysis workbench. 
Bioinformatics 25: 1189-1191. 

60. Cole C, BarberJD, Barton C>] (2008) I'he Jpred 3 secondary structure prediction 
ser\'er. Nucleic aeids res 36: \V197-\V201. 

61. Wheclan S, (joldman N. (2001) A general empirical model of protein evolution 
derived from multiple protein families using a maximum-likelihood approach. 
Mol Biol Evol 18(5): 691-9. 

62. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (201 1) MEGA5: 
molecular evolutionary genetics analysis using maximum likelihood, evolution- 
ary distance, and maximum parsimony methods. Mol Biol Evol 28: 2731-2739. 

63. Mortazavi A, Williams BA, Mecue K, Schacffer L, Wold B (2008) Mapping and 
quantifying mammalian transcriptomes bv RNA-Seq. Nat Methods 5: 1—8. 

64. Guo S, KimJ (2007) Molecular evolution of Drosophila odorant receptor genes. 
Mol Biol Evol 24(5): 1198-207. doi:10.1093/molbev/msm038 

65. Benton R (2006) On the ORigin of smell: odorant receptors in insects. Cell Mol 
Life Sci 63: 1579-1585. 

66. Pitts RJ, Fox AN, Zwiebel LJ (2004) A highly conserved candidate 
ehemoreceptor expressed in both olfactory and gustatory tissues in the malaria 
vector Anopheles gambiae. Proc Nati Acad Sci USA 101: 5058-5063. 

67. Bouyer J, Pruvot M, Bengaly Z, Guerin PM, Lancelot R (2007) Learning 
influences host choice in tsetse. Biology letters 3: 113—116. doi:10.1098/ 
rsbl.2006.0578. 

68. Hargrove JW (1976) Amino acid metabohsm during flight in tsetse flies. J Insect 
Physiol 22(2): 309-313. 

69. Bellmann D, Riehardt A, Freyberger R, Nuwal N, Schwiirze! M, et al. (2010) 
OptogcneUcally Induced Olfactory Stimulation in Drosophila Larvae Reveals 
the Neuronal Basis of Odor-Aversion behavior. Front Behav Neurosci 4(June): 
27. doi:10.3389/fnbeh.2010.00027 

70. Marinotti O, Calvo E, Nguyen QK, Dissanayake S, RibeiroJM, et al. (2006) 
Gcnomc-wide analysis of gene expression in adult Anopheles gambiae. Insect Mol 
Biol 15: 1-12 

71. Parisi M, Nuttall R, Edwards P, Minor J, Naiman D, et al. (2004) A survey of 
ovary-, testis-, and soma-biased gene expression in Drosophila melanogaster adults. 
Genome Biol 5: R40. 

72. McGraw L.\, (?.lark AG, Wolfner MF (2008) Post-mating gene expression 
profiles of female Drosophila melanogaster in response to time and to four male 
accessory gland proteins. Genetics 179: 1395—1408. 

73. Bionaz M, Periasamy K, Rodriguez-Zas SL, Everts RE, Lewin HA, et al. (2012) 
Old and new stories: revelations from functional analysis of the bovine 
mammary transcriptome during the lactation cycle. PLoS ONE 7: e33268. 

74. Bohbot J, Pitts RJ, Kwon HW, Rutzler M,' Robertson HM, et al. (2007) 
Molecular characterization of the Aedes aegj^pti odorant receptor gene family. 
Insect Mol Biol 16: 525-537. 

75. Zhou X, Slone JD, Rokas A, Bcrger SL, LiebigJ, et al (2012) Phylogenetic and 
Transcriptomic Analysis of Chemosensory Receptors in a Pair of Divergent Ant 
Species Reveals Sex-Specific Signatures of Odor Coding. PLOS Genet 8(8): 
el002930. 

76. Hill CA, Fox AN, Pitts RJ, Kent LB, Tan PL,et al (2002) G Protein-Coupled 
Receptors in. Anopheles gambiae. Science 298: 176—178. 

77. Kent LB, Walden KKO, Robertson HM (2008) The Gr Family of Candidate 
Gustatory and Olfactory Receptors in the Yellow-Fever Mosquito Aedes aegypti. 
Chem Senses 33: 79 93. 

78. Richards S, Gibbs R, VVeinstoek C;, Brown S, Danell R et al (2008) The genome 
of the model beetle and pest TrihoUum casteneum. Nature 452(7190): 949-955. 

79. Engsontia P, Sanderson AP, Cobb M, Waldcn KKO, Robertson HM et al 
(2008) The red flour beetie's large nose: An expanded odorjuit receptor gene 
family in Tribolium casteneum. Insect Biochemistry and Molecular Biology 38: 
387-397. 

80. Robertson HM, Gadau J, Wanner KW (2010) The insect chemoreecptor 
superfamily of the parasitoid jewel wasp JVasonia vitripennis. Insect Mol Biol 19(1): 
121-136. 



PLOS Neglected Tropical Diseases | www.plosntds.org 



8 



April 2014 | Volume 8 | Issue 4 | e2663 



